The thermodynamic Casimir effect in the neighbourhood of the λ-transition: A Monte Carlo study of an improved three dimensional lattice model

We study the thermodynamic Casimir effect in thin films in the three dimensional XY universality class. To this end, we simulate the improved two component φ4 model on the simple cubic lattice. We use lattices up to the thickness L0 = 33. Based on the results of our Monte Carlo simulations we compute the universal finite size scaling function θ that characterizes the behaviour of the thermodynamic Casimir force in the neighbourhood of the critical point. We confirm that leading corrections to the universal finite size scaling behaviour due to free boundary conditions can be expressed by an effective thickness L0,eff = L0 + Ls, with Ls = 1.02(7). Our results are compared with experiments on films of 4He near the λ-transition, previous Monte Carlo simulations of the XY model on the simple cubic lattice and fieldtheoretic results. Our result for the finite size scaling function θ is essentially consistent with the experiments on films of 4He and the previous Monte Carlo simulations.